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ABSTRACT 



We investigate the distribution of energy density in a stationary self-reproducing 
inflationary universe. We show that the main fraction of volume of the universe in 
a state with a given density p at any given moment of proper time t is concentrated 
near the centers of deep exponentially wide spherically symmetric wells in the density 
distribution. Since this statement is very surprising and counterintuitive, we perform 
our investigation by three different analytical methods to verify our conclusions, and 
then confirm our analytical results by computer simulations. If one assumes that we 
are typical observers living in the universe at a given moment of time, then our results 
may imply that we should live near the center of a deep and exponentially large void, 
which we will call infloid. Validity of this particular interpretation of our results is 
not quite clear since it depends on the as-yet unsolved problem of measure in quan- 
tum cosmology. Therefore at the moment we would prefer to consider our results 
simply as a demonstration of nontrivial properties of the hypersurface of a given time 
in the fractal self-reproducing universe, without making any far-reaching conclusions 
concerning the structure of our own part of the universe. Still we believe that our 
results may be of some importance since they demonstrate that nonperturbative ef- 
fects in quantum cosmology, at least in principle, may have significant observational 
consequences, including an apparent violation of the Copernican principle. 



1 Introduction 



According to the Copernican principle, the only special thing about the Earth is that we are 
living here. We are not at the center of the universe, as people thought before. This point 
of view is reflected also in the so-called cosmological principle, which asserts that our place 
in the universe is by no means special and that the space around us has to be homogeneous 
and isotropic after smoothing over small lumps of matter. This principle lies in the foundation 
of contemporary cosmology Jl| since it has not only definite philosophical appeal but also an 
apparent observational confirmation by a host of data on large scale structure of the universe. 
However, theoretical interpretation of this principle is usually based on the big bang picture of 
the universe and its evolution, inherently related to simple geometry of Friedmann-Robertson- 
Walker type. The only theoretical justification of homogeneity and isotropy of the universe which 
is known to us at present is based on inflationary cosmology. But this theory simultaneously with 
explaining why our universe locally looks so homogeneous predicts that on an extremely large 
scale the universe must be extremely inhomogeneous @. Thus, after providing certain support 
to the cosmological principle, inflationary theory eventually removes it as having only limited 
validity. But until very recently we did not suspect that inflation may invalidate the Copernican 
principle as well, since there is nothing about inflation which would require us to live in the center 
of the universe. 

The situation became less obvious when we studied the global structure of inflationary universe 
in the chaotic inflation scenario, and found that according to a very wide class of inflationary 
theories, the main fraction of volume of the universe in a state with a given density p at any 
given moment of time t (during or after inflation) should be concentrated near the centers of 
deep exponentially wide spherically symmetric wells in the density distribution ||. This result 
is based on investigation of nonperturbative effects in the theory of a fractal self-reproducing 
universe in the chaotic inflation scenario [|J.[] 

Observational implications of this result depend on its interpretation. If we assume that we 
live in a part which is typical, and by "typical" we mean those parts of the universe which have 
the greatest volume with other parameters (time and density) being equal, then our result implies 
that we should live near the center of one of the wells in the density distribution. There should 
be many such wells in the universe, but each of them should be exponentially wide. In what 
follows we will call these wells "infloids." An observer living near the center of an infloid will see 
himself "in the center of the world," which would obviously contradict the Copernican principle. 

One should clearly distinguish between the validity of our result and the validity of its inter- 
pretation suggested above. Even though the effect by itself is rather surprising we think that it is 
correct. We verified its validity by three independent analytical methods, as well as by computer 
simulations. Meanwhile the validity of its interpretation is much less clear. The main problem is 
related to the ambiguity in the choice of measure in quantum cosmology 0. There are infinitely 

1 Self-reproduction of the universe is possible in the new inflationary theory as well ||, but as we will see, in 
this theory the effect which we are going to discuss is negligibly small. 
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many domains with similar properties in a self-reproducing inflationary universe. When we are 
trying to compare their volumes, we are comparing infinities. The results of this comparison 
depend on the choice of the regularization procedure. The prescription that we should compare 
volumes at a given time t in synchronous coordinates is intuitively appealing, but there exist other 
prescriptions which lead to different conclusions || [7|, ||. Until the interpretation problem is 
resolved, we will be unable to say for sure whether inflationary cosmology actually predicts that 
we should live in a center of a spherically symmetric well. Still this possibility is so interesting 
that it deserves a detailed investigation even at our present, admittedly rather incomplete level 
of understanding of quantum cosmology. This is the main purpose of our paper. 

In Section |2| we will give a short review of the theory of self-reproducing universe in the chaotic 
inflation scenario and discuss which type of phenomena should be called typical in such universe. 
Then we will describe two approaches to the problem of estimating the typical magnitude of 
the quantum fluctuations under the volume weighted measure. The first is based on counting 
the balance of probability factors. The second is based on the investigation of the probability 
distribution P p (4>, t). This distribution describes the portion of the physical volume of the universe 
which contains the field <fi at the time t. According to ||, this distribution rapidly approaches 
a stationary regime, where the portion of the physical volume of the universe containing the 
field <p becomes independent on time. Investigation of this distribution in Section |3| will allow 
us to derive our result in a different way. In Section [5] we will develop a path integral approach 
to the investigation of P p (<p,t). The new method provides another way to confirm our results. 
However, this method is interesting by itself. It gives us a new powerful tool for investigation 
of the global structure of the self-reproducing universe, which may be useful independently of 
existence of the effect discussed in this paper. In Section |6| we will describe computer simulations 
which we used to verify our analytical results. Only then, after we make sure that our rather 
counterintuitive results are actually correct, we will describe their possible interpretation and their 
observational consequences. In Section [7] we will describe the structure of infloids, their evolution 
after the end of inflation and their observational manifestations. In Section || we will discuss 
our results, ambiguities of their interpretation, and formulate our conclusions. In Appendix we 
present generalization of our results for different time parametrizations. 



2 The Self-Reproducing Universe 

Let us consider the simplest model of chaotic inflation based on the theory of a scalar field 
minimally coupled to gravity, with the effective potential V(4>). If the classical field <fi (the inflaton 
field) is sufficiently homogeneous in some domain of the Universe, then its behavior inside this 
domain is governed by the equations 

+ 3#0 = -U'(0) , (1) 
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Here H = a/a, a(t) is the scale factor of the Universe, k = +1, —1, or for a closed, open or flat 
Universe, respectively. M p is the Planck mass, which we will put equal to one in the rest of the 
paper. 

Investigation of these equations has shown that for many potentials V(<f)) (e.g., in all power- 
law V(4>) ~ <p n and exponential V((p) ~ e°"^ potentials) there exists an intermediate asymptotic 
regime of slow rolling of the field <fi and quasi-exponential expansion (inflation) of the Universe 
|2|. At this stage, which is called inflation, one can neglect the term in (|2|), as well as the terms 
Jj- and ^-<fi 2 in @. Therefore during inflation 




H = \l—, *=-^h- 0) 
In the theories V(4>) ~ 4> n inflation ends at <p = <ft e , where <p e ~ 10~ 1 n. 

Inflation stretches all initial inhomogeneities. Therefore, if the evolution of the Universe were 
governed solely by classical equations of motion, we would end up with an extremely smooth 
Universe with no primordial fluctuations to initiate the growth of galaxies. Fortunately, new 
density perturbations are generated during inflation due to quantum effects. The wavelengths of 
all vacuum fluctuations of the scalar field <fi grow exponentially in the expanding Universe. When 
the wavelength of any particular fluctuation becomes greater than if -1 , this fluctuation stops 
oscillating, and its amplitude freezes at some nonzero value 5<p(x) because of the large friction 
term 3H(p in the equation of motion of the field <fi. The amplitude of this fluctuation then remains 
almost unchanged for a very long time, whereas its wavelength grows exponentially. Therefore, 
the appearance of such a frozen fluctuation is equivalent to the appearance of a classical field 
S(j)(x) that does not vanish after averaging over macroscopic intervals of space and time. 

Because the vacuum contains fluctuations of all wavelengths, inflation leads to the creation 
of more and more perturbations of the classical field with wavelengths greater than H . The 
average amplitude of such perturbations generated during a time interval if -1 (in which the 
Universe expands by a factor of e) is given by 



IWI^. (4) 

The phases of each wave are random. It is important also that quantum fluctuations occur 
independently in all domains of inflationary universe of a size greater than the radius of the event 
horizon H~ l . Therefore, the sum of all waves at any given region of a size 0(H~ l ) fluctuates 
and experiences Brownian jumps in all directions in the space of fields. The standard way 
of description of the stochastic behavior of the inflaton field during the slow-rolling stage is to 



coarse-grain it over separate domains of radius H (we will call these domains "/i-regions" R 110 
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to indicate that each of them has the radius coinciding with the radius of the event horizon H : ) 
and consider the effective equation of motion of the long-wavelength field [11. 



d V'(cf>) H*I\<P) 

Here £(t) is the effective white noise generated by quantum fluctuations. 

Let us find the critical value 0* such that for V{<p) < V(0#) the classical slow roll dominates the 
evolution of the inflaton, while for V(4>) > V(0*) the quantum fluctuations are more important. 
Within the characteristic time interval At = H^ 1 for values of inflaton near the critical value 0* 
the classical decrease A0 = At of the inflaton, defined through (||), is of the same magnitude 
as the typical quantum fluctuation generated during the same period, given by (0j). After some 
algebra we get from equations (|3|) and (§) the relation defining 0* implicitly: 

3H 3 (4>.) 41/(0.) . ,„> 



Let us consider for definiteness the theory V((p) = A0 4 /4. In this case equation (f| yields 
0* ~ A -1 / 6 . One can easily see that if < 0*, then the decrease of the field due to its classical 
motion A0 = 1/2-7T0 is much greater than the average amplitude of the quantum fluctuations 

50 = ^/a/67t0 2 generated during the same characteristic time interval H~ x . But for > 0*, 
8(j)(x) will exceed A0, i.e. the Brownian motion of the field will become more rapid than its 
classical motion. Because the typical wavelength of the fluctuations 5<fi(x) generated during this 
time is H~ l , the whole Hubble domain after the time H~ l becomes effectively divided into e 3 
/i-regions, each containing almost homogeneous (but different from each other) field — A0 + 50. 

In almost half of these domains (i.e. in e 3 /2 ~ 10 /i-regions) the field grows by \S<j)(x) \ — A0 w 
|50(a;)| = H/2tt, rather than decreases. During the next time interval At = H^ 1 the field grows 
again in the half of the new /i-regions. Thus, the total number of /i-regions containing growing 
field becomes equal to (e 3 /2) 2 = e 2 ( 3-ln2 ). This means that until the fluctuations of field grow 
sufficiently large, the total physical volume occupied by permanently growing field (i.e. the total 
number of /i-regions containing the growing field 0) increases with time like exp[(3 — In 2) Ht]. 
This leads to the self-reproduction of inflationary domains with > 0* in the chaotic inflation 
scenario 



Note, that the greater is the value of the effective potential, the greater is the rate of ex- 
ponential expansion of the universe. As a result, the main growth of the total volume of the 
universe occurs due to exponential expansion of the domains with the greatest possible values of 
the Hubble constant H = H mSuX 0, H. In some models there is no upper bound to the value of 



H [12, Q3l. However, in the simplest versions of chaotic inflation based on the Einstein theory 



of gravity there are several reasons to expect that there exists an upper bound for the rate of 
inflation @, H [Hj . 
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In what follows we will assume that there is an upper bound H max on the value of the Hubble 
constant during inflation. For definiteness we will assume that H max = yW/3, which corresponds 
to the Planck boundary V(<p p ) = 1. This is a rather natural assumption for chaotic inflation. 
However, one should note that in some models H max may be much smaller. In particular, in the 
new inflation scenario H max = ^j8nV(0)/3 is many orders of magnitude smaller than 1. 

The independence of the subsequent evolution of the /i-region on its previous history, the 
dominance of the domains where the inflaton field energy grows rather than decreases in the 
volume weighted measure and the upper bound for the energies at which the inflation can proceed 
are the three main features inherent to many models of inflation. When all these features are 
present the evolution of the inflationary universe as a whole approaches regime which we called 
global stationarity in JJ. This stage is characterized by the stability of the distribution of regions 
with various local values of energy density and other parameters, while the number of such 
regions grows exponentially with a constant coefficient, proportional to the maximal possible 
rate of inflation Ai = df T ff mH . Here is a model dependent fractal dimension of the classical 
space [||, P , which is very close to 3 for small coupling constants of the inflaton field. 

The new picture of the universe is extremely unusual, and it may force us to reconsider our 
definition of what is typical and what is not. In particular, the standard theory of the large scale 
structure of the universe is based on the assumption that a typical behavior of the scalar field at 
the last stages of inflation is described by equations (||), (|). This is indeed the case if one studies 
a single branch of inflationary universe beginning at <C 0*. However, if one investigates the 
global structure of the universes at all and tries to find the typical behavior of all inflationary 
domains with a volume weighted measure, the result may appear to be somewhat different. 



3 Stationary Inflation and Nonperturbative Effects 

Suppose that we have one inflationary domain of initial size H~ 1 , containing scalar field > 0*. 
Let us wait 15 billion years (in synchronous time t in each part of this domain) and see what are 
the typical properties of those parts of our original domain which at the present moment have 
some particular value of density, e.g. p = 10~ 29 g ■ cm -3 . The answer to this question proves to 
be rather unexpected. 

This domain exponentially expands, and becomes divided into many new domains of size 
iy -1 , which evolve independently of each other. In many new domains the scalar field decreases 
because of classical rolling and quantum fluctuations. The rate of expansion of these domains 
rapidly decreases, and they give a relatively small contribution to the total volume of those parts 
of the universe which will have density 10~ 29 g • cm -3 15 billion years later. Meanwhile those 
domains where quantum jumps occur in the direction of growth of the field gradually push this 
field towards the upper bound where inflation can possibly exist, which is presumably close to 
the Planck boundary V(0 P ) ~ 1. Such domains for a long time stay near the Planck boundary, 
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and exponentially grow with the Planckian speed. Thus, the longer they stay near the Planck 
boundary, the greater contribution to the volume of the universe they give. 

However, the domains of interest for us eventually should roll down and evolve into the 
regions with density 10~ 29 g ■ cm -3 . Thus, these domains cannot stay near the Planck boundary 
for indefinitely long time, producing new volume with the Planckian speed. However, they will 
do their best if they stay there as long as it is possible, so that to roll down at the latest possible 
moment. In fact they will do even better if they stay near the Planck boundary even longer, to 
save time for additional rapid inflation, and then rush down with the speed exceeding the speed 
of classical rolling. This may happen if quantum fluctuations coherently add up to large quantum 
jumps towards small 0. This process is dual to the process of perpetual climbing up, which leads 
to the self-reproduction of inflationary universe. 

Of course, the probability of large quantum jumps down is exponentially suppressed. However, 
by staying longer near the Planck boundary inflationary domains get an additional exponentially 
large contribution to their volume. These two exponential factors compete with each other to give 
us an optimal trajectory by which the scalar field rushes down in those domains which eventually 
give the leading contribution to the volume of the universe. From what we are saying it should 
be clear that the quantum jumps of the scalar field along such optimal trajectories should have 
a greater amplitude than their regular value J^, and they should preferably occur downwards. 
As a result, the energy density along these optimal trajectories will be smaller than the energy 
density of their lazy neighbors which prefer to slide down without too much of jumping. This 
creates wells in the distribution of energy density, which we called infloids J3|. 

Suppose that the extra time interval spent at highest energies is At. Then we win the volume 
by factor of exp (d{ v H max Atj . However, to compensate for the lost time the inflaton field has to 

jump at least once (let us say, when it reaches the value 0) with the amplitude 50 = n(0) if (0)/27r 
such that it covers in one jump the distance which would otherwise require time At to slowly roll 
down: 



& ( ,)4 = !« = bW «™ (7) 

V W) 

where we introduced the factor n(0) by which the jump is amplified, i.e. by which it is greater 
than the standard jump if(0)/27r. The probability of such jump is suppressed by the factor 
exp |n 2 (0)). The leading contribution to the volume of the universe occurs due to the jumps 
which maximize the volume weighted probability: 



P ~ exp \d b H m3X At(4>) - -n (0)J = exp \ d b H maK n{4>) ~ g 71 (0) J • ( 8 ) 

Maximizing with respect to n(0) gives the amplification factor as a function of the location of 
the jump on the inflaton trajectory: 
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n{4>) = 4 dft if max . (9) 

In fact, we have found 0] that the typical trajectories which give the leading contribution to the 
volume of the universe consist entirely of such subsequent jumps. In what follows we will give an 
alternative derivation of this result. Meanwhile, comparing with (^) one immediately sees that 
n(4>) S> 1 for < 0*, since d{ r ~ 3 and H max ^> H((f>) for such values of inflaton field in chaotic 
inflation. Therefore, our treatment of these quantum fluctuations as large and rear quantum 
jumps is self-consistent. 

To avoid misunderstandings one should note that a more accurate definition of amplification 
coefficient would be n((j)) + 1. Indeed, in the absence of nonperturbative effects we would have 
n((j)) = since perturbative jumps occur in both directions with equal probability. The coefficient 
n((f>) relates an additional amplitude of jumps down to the regular perturbative amplitude of the 
jumps in both directions. This subtlety will not be important for us here since we are interested 
in the case n ^> 1. 

It is interesting that the coefficient of amplification n(<f>) can be directly related to the ratio of 
amplitudes of conventional scalar and tensor perturbations generated at the same scale at which 
the jump occurs. The amplitudes of these perturbations can be written as follows: 

4rw=(!e) -c s m, Arw = ( s A =e T ^. do) 

\P/s 2vr0 \PJt M p 

Here c$ and ct are some coefficients of the order of unity. Using these expressions we can rewrite 
(§) for G?f r ~ 3 in the form: 



3c T H max A p s ert (4>) 

In the same way as the conventional amplitude of jumps H/2n is related with the well known 
perturbations of the background energy density, the "nonperturbatively amplified" jumps which 
we have just described are related to the "nonperturbative" contribution to deviations of the 
background energy density from its average value. A possible interpretation of this result is that 
at the length scale associated with the value of the field <fi there is an additional nonperturbative 
contribution to the monopole amplitude: 



„^ ) = g^|gg^ r . (12) 

We will discuss the structure of infloids and their possible observational consequences in 
Section |7|. Here we only note that eq. ( |TT| ) gives a simple tool for understanding of the possible 
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significance of the effect under consideration. Indeed, in the simplest chaotic inflation models, 
such as the theory -</> n , one has H max ~ M p and A p s ert ((j)) ^> A^p Tt ((f)); thus one has n(<f>) ^> 1. 
On the other hand, in the versions of chaotic inflation scenario where inflation occurs near a 
local maximum of the effective potential (as in the new inflation models) -ff max is many orders 
of magnitude smaller than M p , and therefore the non-perturbative effects discussed above are 
negligibly small. Thus, investigation of nonperturbative effects can give us a rather unexpected 
possibility to distinguish between various classes of inflationary models. We will return to this 
issue in the end of the paper. 



4 Nonperturbative effects and branching diffusion 



One of the best ways to examine nonperturbative effects is to investigate the probability distri- 
bution P p ((j), t) to find a domain of a given physical volume in a state with a given field at some 
moment of time t. The distribution P p ((f),t) obeys the following branching diffusion equation 
PL H6|, IBIl: 



dt d<f) \ 2V2tt <90 V 2y/2n P J 3H(</>) P J KV> 7 1 J 

This equation is valid only during inflation, which typically occurs within some limited interval 
of values of the field <fi: <f) m in < <p < 4> m ax- in the simplest versions of chaotic inflation model 
<j) m i n = <p e ~ 1, where <p e is the boundary at which inflation ends. Meanwhile, as we argued in the 
previous section, 4> ma x is close to the Planck boundary <p p , where V((f> p ) = 1. To find solutions of 
this equation one must specify boundary conditions. Behavior of solutions typically is not very 
sensitive to the boundary conditions at e ; it is sufficient to assume that the diffusion coefficient 
(and, correspondingly, the double derivative term in the r.h.s. of equation flilf )) vanishes for 
< e M. The conditions near the Planck boundary play a more important role. In this paper 
we will assume, that there can be no inflation at V(<j>) > 1, which corresponds to the boundary 
condition P p (4>, t)\<j>><t> p = 0. In the end of the paper we will discuss possible modifications of our 
results if 4> max differs from cf) p . 

One may try to obtain solutions of equation ( |T5D in the form of the eigenfunction series: 

oo 

P P (<M) = £^(0) e^n^), (14) 

s=l 

where, in the limit of large time t, only the term with the largest eigenvalue Ai survives. The 
function 7Ti(0) in the limit t — ► oo has a meaning of a normalized time-independent probability 
distribution (so called invariant probability density of the branching diffusion) to find a given 
field in a unit physical volume, whereas the factor e Xlt shows the overall growth of the volume 
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of all parts of the universe, which does not depend on <p in the limit t — > oo. This "ground state" 
eigenfunction satisfies the following equation: 



In the limit when we can neglect the diffusion (second derivative) term it is easy to solve this 
equation: 



vti(0) = C(0 O ) y,^ exp 



00 r 



Ai 



3ff(Q 9ff 2 (C) 



(16) 



where we chose some starting point O and the corresponding normalization constant C(4>o) 
which should match this approximate solution to the exact one at this point. As before, let 
us introduce the fractal dimension of classical space-time through Ai = df r H max (see || H for 
detailed discussion of the fractal structure of self- reproducing universe). Let us also introduce 
the critical value of inflaton 0f r at which the no-diffusion approximation for flT5|) breaks. Then, 
since H maiX ^> H((j)) for chaotic inflation, we can rewrite (IB) as: 



tti(0) = C(0fr) v exp 



V 



(17) 



Substituting (17) into fll5|) we get the defining relation for the value of inflaton field 0f r at 
which the no-diffusion approximation breaks: 



Ai 



9^ 5 (0fr) 

4tt 2 {V'{4>b)f 



1 . 



;is) 



We can rewrite (the square root of) this relation in a form which makes the comparison with 
the definition of the other critical value 0* more apparent: 



V #(</> fr ) 27rF'(0 fr ) ' 1 ; 

Comparing (|19D with (|6]), one finds that for all chaotic inflation models 0f r < 0* (one can as- 
sume self-consistently that H max ^> H(<p{ r ) in such models). The value of 0* in (D) comes from 
comparing the slow roll rate in a given /i-region with the typical amplitude of quantum fluctua- 
tions while considering only the /i-regions generated locally from the region which we picked. On 
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the other hand, the value 0f r comes from comparing the slow roll rate to the typical amplitude 
of fluctuations considering all /i-regions in the whole universe which happen to have the same 
value of inflaton field inside. The fact that the second constraint is more stringent is yet another 
indication of the considerably larger magnitude of the quantum fluctuations when we take into 
account the whole stationary universe. 

In the particular case of the simplest theory with V(0) = A0 4 /4, we have H = yj2n\/3(f) 2 , 
0fr ~ A -1 / 8 <C 0*, and the dependence of the solution (|17|) on is [§]: 



7n(0) ~0Vx^ ff — . (20) 

This is an extremely strong dependence. For example, for the realistic value of the coupling 
constant A ~ 1CT 13 chosen to fit the observable large scale structure of the universe one has 
cff r ~ 3. One may assume for definiteness that H max = y^87r/3, corresponding to inflation with 

V{4>) = 1 (Planck density). Then one has an extremely sharp dependence n\ ~ 1 - 2 ' 108 . All 
surprising results we are going to obtain are rooted in this effect. One of the consequences is the 
distribution of energy density p. For example, during inflation p « A0 4 /4. Equation (|20|) implies 
that the distribution of density p is 

Pp{p) ~ P 3 ' 107 • (21) 

Thus at each moment of time t the universe consists of indefinitely large number of domains 
containing matter with all possible values of density, the total volume of all domains with density 
2p being approximately 10 10 times greater than the total volume of all domains with density p\ 

Let us consider now all inflationary domains which contain a given field at a given moment 
of time t. One may ask the question, what was the value of this field in those domains at the 
moment t — if -1 ? In order to answer this question one should add to the value of its classical 
drift (pH^ 1 and the amplitude of quantum jumps A<f). The typical jump is given by <50 = ±H/2ir. 
At the last stages of inflation this quantity is by many orders of magnitude smaller than 0i/~ 1 . 
But in which sense jumps ±H/2tt are typical? If we consider any particular initial value of the 
field 0, then the typical jump from this point is indeed given by ±H/2n under the conventional 
comoving measure. However, if we are considering all domains with a given and trying to find 
all those domains from which the field could originate back in time, the answer may be quite 
different. Indeed, the total volume of all domains with a given field at any moment of time t 
depends on extremely strongly: the dependence is exponential in general case (|TTD, or a power 
law with a huge power, like in the case of A0 4 /4 theory (^0|). This means that the total volume of 
all domains which could jump towards the given field from the value + A0 will be enhanced 
by a large additional factor P p (0 + A0)/P p (0). On the other hand, the probability of large jumps 
A0 is suppressed by the Gaussian factor exp (— 2n 2 (A0) 2 /if 2 ). Thus, under the established 
stationary probability distribution the probability of the inflaton field in a given domain to have 
experienced a quantum jump A0 is given by: 
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WHeip U„„^.^>!), (22) 



One can easily verify that this distribution has a sharp maximum at: 



A*, = *«-^^=«(«^, (23) 

and the width of this maximum is of the order In other words, most of the domains of a given 
field are formed due to non-perturbative (hence the subscript "np") jumps which are greater 
than the "typical" ones by a factor n(4>) which coincides with our previous result @. For future 
reference, we will write here this result in an equivalent form, 

= < 24 > 

The limit of applicability of this expression is below the energy level V(0f r ) (see fll8|), (|19"1) for 
definition of the critical value 0f r ). 

In particular, for the theory A0 4 /4 we have 

n(0) = Ai0 . (25) 



For iJ max = y 87r/3, A 1 and ~ 4.5, which corresponds to today's horizon scale, this gives 
the amplification coefficient 

n{(j)) = 2 V / 6V0 ~ 40 . (26) 



5 Volume Weighted Slow Rolling Approximation 



We learned in the previous section that quantum fluctuations in volume weighted measure have 
pretty large expectation value, which makes the jumps to go preferentially downwards (unlike in 
comoving measure where there is no preferred direction of the fluctuations and therefore they 
have zero expectation value). Indeed, such was the very meaning of our derivation of large jumps 
that they had to occur in the direction of usual slow roll in order to make up the extra time 
spent by inflaton at higher energies. Therefore, we can conclude that the slow rolling speed itself 
gets a correction corresponding to the rate at which such large jumps occur and their size. Since 
each such jump of the size n(4>) ■ H(<fi)/2it occurs during time interval if~ 1 (0), we can estimate 
the additional speed gained by the inflaton as n(<f>) ■ H 2 ((j))/2ir, thus bringing the overall slow roll 
speed to the volume weighted value (we substituted @ for the value of n((f>), the amplification 
factor): 
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dfr ffmax ou ,/^ • (27) 



3 #(</>) Ir mdX 3V'(<p) 

Here the minus sign in front of the correction term is due to the preferred direction of the jumps 
bringing the slow roll speed to higher absolute value. 



The limits of applicability of this expression are the same as for fl23|), i. e. below energy density 
corresponding to the critical value 0f r of inflaton field, defined by (|18|), ([TTj|). However, those limits 
simply tell where the approximate expression fl27|) is valid, while the effect of speeding up the 
slow roll of the inflaton is valid in a much wider range. 

Let us derive a more general version of this result and, correspondingly, a more general 



expression for amplified quantum jumps @, (|23|) which will be valid for almost whole range of 



variation of the inflaton field. The volume weighted probability distribution can be defined as 
the path integral over all realizations of noise taken with gaussian weight modified by the volume 
factor HTn, HI: 



P p {<f>,t) =Jv£expU (-^ 2 (s)+3#(<^(s))) ds j HMt)-<f>) ■ (28) 

Here 4>^{s) is the solution of @ with a particular realization of the noise. The gaussian path 
integral over the noise can be converted into the path integral over the histories of inflaton 
evolution if we express the noise through concurrent value of inflaton <ft(t) using the equation 
of motion (Hj: 



2vr • 27TV'((b) . . 



It is convenient to make the following variable transformation: 



277 - >' . (30) 



# 3 / 2 (0') ^ ' 

In terms of this variable the definition of the white noise is rewritten in compact form 



C(t) = -z + W{z) , (31) 

where we introduced the "superpotentiafQ (we used the relation (|3H| ) to re-express it in terms of 
the derivative with respect to z): 

2 This name is due to the fact that W(z) plays a role of a superpotential in a SUSY-Schrodinger like version of 
Fokker-Planck equation. 
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2tt d V((j)) d f 3 \ 

W{Z} ~ 3^/2(0) ~ Tz \l6V(zj) ■ (32) 

The path integral defining the volume weighted measure in terms of z(t) becomes, after 
substituting (0) into (|28]): 



P p (z,t) = JVz(s)J[z]exp<- J (^[z(s)~W(z(s))] 2 -3H(z(s))^j ds\ 6 (z(t) - z) . (33) 



The Jacobian J[z] of the transformation from £ to and then to z is pre-exponential jT7 and 
unimportant for our current investigation. We will neglect it in what follows. 

Let us find the trajectory z(t) [which we will translate later into trajectory (f>(t)\ which con- 
tributes most to the path integral (|33|). Such saddle point trajectory will correspond to the typical 
history of evolution of inflaton under volume weighted measure. The exponent in the path integral 
(|33|) looks like a Euclidean version of Lagrangian action, which corresponds to interpretation of 
diffusion equation (|T3|) as a Euclidean Schrodinger equation for a point particle. We can rewrite 
this action in Hamiltonian form using the conventional relation: 



z(t) 



Jcdt= J pdz- Jndt, (34) 



where the canonical momentum is 



dL 

p = — = z- W(z) . (35) 
oz 



Since the action does not contain explicit time dependence, the Hamiltonian is conserved: 

H = ip 2 +pW(z) + 3H(z) = Ai . (36) 

The reason why the conserved Hamiltonian is equal to the highest eigenvalue is that at the end 
we should get the time dependence of a type exp (Ai t) as warranted by the stationary solution 
([14]) . Meanwhile, the J pdz term of the action should give us the correct (semi-classical) field 
dependence of the probability density P p (z((j)),t) (see below). 

Solving the Hamiltonian constraint (|36|) with respect to p (we have to choose the positive 
solution of the equation for rolling down), and using (3"5|), we obtain the equation for the typical 
volume weighted trajectory: 

z = ^JW 2 (z) + 2Ai - 6H(z) . (37) 
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This equation translates back in terms of inflaton field variable into volume weighted slow roll 
equation: 



\ Km*)) 



+ ( dfr H n 



3H{<j>) 



H 3 ((f>) 
2tt 2 



(38) 



For most of the inflaton range of variation (except very close to Planck boundary) we can ignore 
3H term with respect to df r H mSuX term. The relative importance of the two remaining terms 
under the square root is governed by the critical value 0f r - - below this level the first term 
is more important, while above it the second one dominates. Not surprisingly, below 0f r this 
equation coincides with (p7|). However, its validity limits are much wider, allowing us to use it 
beyond 0f r , to which the applicability of ( pT|) was limited. 



We can write down a good approximation for the field dependent normalized probability 
density 7Ti (</>), omitting the less important pre-exponential terms: 



7Tl(0) 



exp < — / pdz 



t 



exp 



V 



3 no ' 



dfr ^max 



3H{<f>) 



8tt 2 



3 HO 

H 3 ((f)) 16V 2 (0 



(39) 



Of course, below 0f r this expression also coincides with its counterpart (|T7|) derived previously. 
This result has remarkable properties which will be studied further in |T8| . 



Using the volume weighted slow roll equation (|38|) we can derive a general expression for the 
amplified quantum jump size. It is given by the change in <p within time interval H~ x calculated 
according to ( pB]) , less the regular expression for the change of field due to the slow roll in comoving 
measure: 



np 



^3^ 2 (0), 



+ (dfr H D 



3H(cf>) 



H(<P) V'{<t>) 



2tt 2 3F 2 (0) 



(40) 



This gives the following expression for the amplification factor (the ration of A0 np and the 
conventional amplitude H/2tt): 



n{</>) 



' 27rV^y 



+ 



2d {r H max - QH{4>) _ 2ttV^) 
H(</>) " 3H^) 



(41) 



The consistency conditions for our results (|38[) - (|4lD arise from several assumptions which 
we made in their derivation and whose validity should be maintained. The first one is that the 
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slow rolling approximation is valid, i.e. <C 3if (0)0. The second is that the amplification factor 
is greater than one. The third condition is that the saddle-point approximation used to derive 
these results is valid, which means that p > 1. And the final, fourth condition is the implicit 
assumption that large quantum jumps which occur in a single /i-region do not make the gradient 
energy inside that region greater than the potential energy of the inflaton field (which, of course, 
would immediately invalidate the inflationary approximation). One can easily check that all four 
conditions lead to the same, very relaxed restrictions — the energy density of the inflaton field 
V{4>) must be lower than the Planck density (or, more precisely, lower than the energy density 
corresponding to the maximal rate of expansion H mauX ). Thus, we can use the results obtained 
above in most of the variation range for the inflaton field in chaotic inflation. 



One can easily check that for < 0f r , H{(f)) <C H max eq. (|4T| ) yields 

n(0) = 4d fr # max ^|| . (42) 

This expression coincides with the expression for the amplification coefficient which we obtained 
earlier by two other methods, see eqs. (|9p and (p3|). 



6 Numerical Simulations 



6.1 The basic idea of computer simulations 

Even though we verified our results by several different methods, they are still very unusual 
and counterintuitive. Therefore we performed a computer simulation of stochastic processes 
in inflationary universe, which allows to obtain an additional verification of our results and to 
calculate the amplification factor n(0) numerically. We have used two different methods of 
computer simulations. The first one is more direct and easy to understand. Its basic idea can be 
explained as follows. 

We have studied a set of domains of initial size H~ l filled with large homogeneous field 0. We 
considered large initial values of 0, which leads to the self- reproduction of inflationary domains. 
From the point of view of stochastic processes which we study, each domain can be modelled by 
a single point with the field in it. Our purpose was to study the typical amplitude of quantum 
jumps of the scalar field in those domains which reached some value 0o = 0(1) close to the end 
of inflation. Then we calculate the amplification factor n(0o) for various 0o- 

Each step of our calculations corresponds to a time change At = uHq 1 . Here H = i7(0 o ), 
and u is some number, u < 1. The results do not depend on u if it is small enough. The evolution 
of the field in each domain consists of several independent parts. First of all, the field evolves 
according to classical equations of motion during inflation, which means that it decreases by "K, 
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during each time interval uHq . Secondly, it makes quantum jumps by 50 = -^yj^fi • Here 
is a set of normal random numbers, which are different for each inflationary domain. To account 
for the growth of physical volume of each domain we used the following procedure. We followed 
each domain until its radius grows two times, and after that we considered it as 8 independent 
domains. In accordance with our condition P p (4>, t)\<j>><f> p = 0, we removed all domains where the 
field jumped to the super-Planckian densities V(<p) > 1. Therefore our method removes the 
overall growth factor e Xlt in the expression P p ~ e Al *7r 1 (0) and directly gives the time-independent 
function 7Ti(0) which we are looking for. Indeed we have checked that after a sufficiently large 
time t the distribution of domains followed by the computer with a good accuracy approached 
the stationary distribution 7i"i (0) which we have obtained in || by a completely different method, 
see Fig. [I]. We used it as a consistency check for our calculations. In what follows we will not 
distinguish between P p and the time- independent factor 7Ti(0). 

We kept in the computer memory information about all jumps of each domain during the 
last time interval Hq 1 before the field inside this domain becomes smaller than 0o- This 
made it possible to evaluate an average sum of all jumps of those domains in which the scalar 
field became smaller than O within the last time interval Hq . Naively, one could expect this 
value to be smaller than since the average amplitude of the jumps is Ssl, but they occur 
both in the positive and negative directions. However, our simulations confirmed our analytical 
result A0 = Ai0 • Tp. In other words, we have found that most of the domains which reach the 
hypersurface = 0o within a time interval At = Hq 1 do it by rolling accompanied by persistent 
jumps down, which have a combined amplitude Ai0o times greater than 

6.2 Details of the method 

Even though this method of calculations may seem quite straightforward, (it is the so-called 
event-tracing Monte-Carlo method) in reality it must be somewhat modified. The main problem 
is obvious if one recalls our expression for the probability distribution P p ~ e Al *7Ti(0) at small 

0: P p ~ e Al * 0V^r Al . As we already mentioned, (omitting the time- dependent factor) this yields 
P p ~ 108 for the realistic value A ~ 1CT 13 . It is extremely difficult to work with distributions 
which are so sharp. 

Therefore in our computer simulations we have studied models with A ~ 0.1, which makes 
computations possible. On the other hand, when one increases the value of A an additional 
problem arises. Our simple expression P p ~ 0v^x Al has been obtained in the limit of very small 
A, which is not perfectly accurate for A ~ 0.1. Therefore we will represent P p in a more general 
form P(0) = (p g ^\ where g(0) approaches a constant value Ai for <C A -1 / 8 . One should 
also take into account that the classical decrease c(0) = of the field during the time 

Hq 1 and the standard deviation s(0) = -^Jj^ (the average amplitude of quantum fluctuations 
during the time Hq 1 ) are not constant throughout the region where the effect takes place. In 
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such a situation an expression for the amplification coefficient n(0) will be slightly different from 
our simple expression n = Ai0. Therefore we should first derive here a more accurate expression 
for n(0), and then compare it with the results of our simulations. 



Consider a point for which = 0o at some time t, at which the stationary regime is already 
established. At the earlier time t — Hq 1 this point was approximately at = O + c(0o) + x, 
where x = 5(p is the sum of all quantum jumps experienced by the field a t this point during 
the last time interval Hq . Consider the probability -P(50) that the field jumped to 0o from 
the point 0o + c(0o) + x - This probability distribution is equal to the distribution P p ((p) times 
the probability of undergoing a quantum fluctuation of length 5(f). 



x 2 



P(x) cx P P (0 O + c(0 o ) + x) exp -— . (43) 



v 2s 2 

Position of the maximum of the distribution P(x) is given by 

P;(0o + c(0 o ) + x) = X_ 
P p (0o + c(0 o ) +x) ~ s 2 ' 



(44) 



To solve this equation for x we need to know P p (0). As earlier we assume P p ((j)) = (f,9 (<!>) w here 
g((f)) varies slowly with cf). If g'<j)hi<j) lng (which happens to be a good approximation for 

<f) ~ 1), eq. ( f44"D can be easily solved, and the expression for n(</>o) looks as follows: 



n ^ = ^T^ ^wTT (^(<Po + c(<p )) 2 + H^)s 2 (cPo) ~ (0o + c(0 o ))) • (45) 
s(0 o ) 2s(0 o ) V v / 

One can obtain a slightly more accurate expression by taking into account dependence of g, c and 
s on 0. Note that in the situation which we are going to investigate O ^ c (0o) ^ x ^ s (0o)- hi 
the limit when g, s, and c can be considered constant, and 4>q ^> s,c this equation leads to our 
earlier expression n(0 o ) = Almo- 
in order to use (|4"5| ) we also need to know g(4>) for our problem. We approximate g(<p) by a 
second order polynomial in and substitute -P(0) = (f ) a o +a ^ +a ^ ( t> into the differential equation for 
P(0). Local analysis around = 0o shows P(0) ~ 56 ~ 23 *~ 4 * . This approximation is accurate 
for ~ 1. 



6.3 Numerical Calculation of n((j)) 

Even for not very small A the distribution P p remains extremely sharp. We have made our 
simulations with A = 0.1, in which case P p ~ 60 . This means that if we want to follow evolution 
of a single domain with = 0.5, then we should simultaneously keep track of 2 60 ~ 10 18 domains 
with = 1. Therefore the simple event-tracing Monte-Carlo approach which we described above 
can be quite adequate for the investigation of P p near its maximum, but not for the study of P p 
far away from the maximum of the distribution. 
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A more advanced approach is to represent the distribution by evenly spaced points with weight 
proportional to the distribution. In other words, we rewrite the probability distribution as a finite 
sum of nearly delta-functional distributions, 



N 

p P {M~Hpi(M > (46) 

i=l 

where 

p .^ t ) = f P P (4>i) for0 l <0<0 l+1 (4?) 
1 otherwise 

At each step of the simulation we investigate the evolution of the distributions pi during the 
time At = uHq 1 . The following equation takes into account classical decrease of the field 0, 
quantum fluctuations, and inflation. 

^^^M-^^hiW) • (48) 

We find P p (4>, t + At) by computing the sum of Pi(<f), t + At). Then we normalize the distribution 
P p ((j),t + At) and again subdivide into a new set of p i: in accordance with (HBp. We repeat this 
process until the resulting distribution P p approaches a stationary regime. 

The most tricky part of the algorithm is to find the amplification factor n(4>o). To do that, 
we associate another distribution Xi(<f),t) with each 0,. Here Xi(<p,t) is the sum of quantum 
fluctuations during the last time interval H , along all trajectories which ended up in the interval 
<pi < <p < <pi+i at the time t. We combine all Xi(<ft,t) into a single distribution X(<f),t), and evolve 
it in the same way as P p (<fi,i), dividing it into the nearly delta-functional distributions Xi(<p,t) 
at every iteration. This is possible because Xi(<f),t) is approximately gaussian and its standard 
deviation is small compared with its mean. When the P p (<j),t) converges, Xo(<f)o,t) approximates 
x, from which n(4>o) can be calculated. 

Decreasing step size u increases accuracy of P p (4>) until some point, after which the accuracy 
starts to decrease. This decrease is explained by the fact that evolved p^s are too sharp and 
therefore are represented inaccurately. To avoid this, having fixed N, we must keep u high 
enough, so that the smallest quantum fluctuation s is wider than the grid spacing. Grid spacing 
is proportional to 1/N and s is proportional to y/u, so the minimal y/u is proportional to 1/N. 
Execution time until convergence is proportional to N 2 / y/u, or for the minimal u it is proportional 
to iV 3 . In practice the largest iV for which the algorithm converges in a reasonable amount of 
time is of the order of 10 3 . 



6.4 Results of Numerical Calculations 



The first step is to verify the numerical algorithm by comparing the probability distribution P(4>) 
it computes with a solution obtained by solving equation (p~5|) obtained in |J. Figure [I] shows 
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Figure 1: Probability distribution P(4>) for V = A0 4 /4, A = 0.1. The dashed line is the numerical 
solution to a differential equation describing P(4>). The solid curve is obtained using computer 
simulations described in this paper. A small deviation between the solid curve and the dashed 
line is due to the finite size of each step and the finite grid size. 

that P{4>) is very close to the correct probability distribution. The deviation between the two 
decreases with step-size. 

The second step is to verify out formula for n. Figure |2] shows that numerically computed 
values of n for different 4>o are close to the ones predicted by the analytical result. The deviation 
is explained by the approximations made in the analytical solution (constancy of g,c,s during 
time H^ 1 ). We have found also, that the typical deviation of the amplitude of jumps from their 
average value nH/2ir is of the order of H/2tt, as suggested by eq. (|22|). This will be important 
for our subsequent considerations. 



7 The Spatial Structure of Infloids 



As one can see from eq. (|38|) , the value of the field <f>{t) corresponding to the typical volume 
weighted trajectories moves down more rapidly that one would expect from the classical slow 



3H(4>) 

trajectories, being surrounded by usual classical neighbors, should correspond to the minima 



roll equation = — 3H ^j\ ■ This is exactly the reason why such nonperturbatively enhanced 
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3.5 - 



2.5 




Figure 2: Comparison between the analytical expression for n(0) (dashed line) and the values 
for n(0) obtained by computer simulations. While the analytical expression is not absolutely 
precise due to various assumptions (such as constancy of g, c, s during the time H' 1 ), it does give 
approximately correct values for n. 



in the distribution of density. To analyse the spatial structure of the universe near the points 
corresponding to the optimal volume-weighted trajectory (38) one should remember that in terms 
of the ordinary comoving measure P c the probability of large fluctuations is suppressed by the 
factor exp (— n 2 (0)/2). It is well known that exponentially suppressed perturbations typically 
give rise to spherically symmetric bubbles \1M . Let us show first of all that the main part of the 
volume of the universe in a state with a given (or with a given density p) corresponds to the 
centers of these bubbles, which we called infloids. 



Consider again the collection of all parts of the universe with a given (or a given density) at 
a given time t. We have found that most of the jumps producing this field during the previous 
time interval if -1 occurred from domains containing the field in a narrow interval of values 
near — <p/H + n(0) • H/2ir. The width of this interval was found to be of the order of H/2tt, 
which is much smaller than the typical depth of our bubble A0 ~ n(<f)) • H/2tt, since we have 
n(0) 3> 1 for all chaotic inflation models. Now suppose that the domain containing the field 
appears not at the center of the bubble, but at its wall. This would mean that the field near the 
center of the bubble is even smaller than 0. Such a configuration could be created by a jump 
from — (j)/H + n(<f)) ■ H/2tt only if the amplitude of the jump is greater than n(0) • H/2tx. 
However, we have found that the main contribution to the volume of domains with a given 
is produced by jumps of an amplitude (n(0) ± 1) • H/2k, the greater deviation from the typical 
amplitude n(0) ■ H/2ir being exponentially suppressed. This means that the scalar field can 
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differ from its value at the center of the bubble by no more than the usual amplitude of scalar 
field perturbations H/2ir, which is smaller than the depth of the bubble by a factor Thus, 
the main fraction of the volume of the universe with a given (or with a given density of matter) 
can be only slightly outside the center. This may lead to a small contribution to anisotropy of 
the microwave background radiation. 

We should emphasize that all our results are based on the investigation of the global structure 
of the universe rather than of the structure of each particular bubble. This is why we assert that 
our effect is non-perturbative. If one neglects that the universe is a fractal and looks only at one 
particular bubble (i.e. at the one in which we live now), then one can find that inside each bubble 
there is a plenty of space far away from its center. Therefore one could conclude that there is 
nothing special about the centers of the bubbles. However, when determining the fraction of 
domains near the centers we were comparing the volumes of all regions of equal density at equal 
time. Meanwhile, the density p wa n of matter on the walls of a bubble is greater than the density 
p cen t er in its center. As we have emphasized in the discussion after eq. (p0|) , the total volume of 
all domains of density p wa n is greater than the total volume of all domains of density p ce ntcr by 
the factor (p W aii/Pcenter) 3 ' 10 • Thus, it is correct that the volume of space outside the center of the 
bubble is not smaller than the space near the center. However, going outside the center brings 
us to the region of a different density, p wa u > p cent er- Our results imply that one can find much 
more space with p = p wa n not at the walls of our bubble, but near the centers of other bubbles. 

This situation can be very schematically illustrated by Fig. [| We do not make an attempt 
to show the spatial distribution of infloids. Rather we show the density distribution near the 
center of each of them. All these regions basically are very similar, but at any particular moment 
of time t there are much more regions with large density since they appeared from the regions 
which inflated at the nearly Planckian density for a longer time. With time the whole set of curves 
should go lower, to smaller p. However, at each moment of time there will be domains with all 
possible values of p, so that the distribution of all curves does not change in time (stationarity). 
If one looks at the whole picture without discriminating between states with different values of 
density, it may seem that there is much more space outside of the centers of the bubbles. However, 
at any given moment of time t the main fraction of volume of the universe in a state with a given 
density p is concentrated near the centers of spherically symmetric bubbles. One may look, for 
example, at the density corresponding to the centers of the third row of curves. At this density 
one may live either near the center of any of the eleven infloids, or at the walls of only three of 
them. The fraction of the volume near the centers would be much greater if we try to show the 
realistic distribution P p (p) ~ p 3 ' 10 of the number of domains with a given density in the theory 
A0 4 /4. 

The nonperturbative jumps down should occur on all scales independently. One may visualize 
the whole process as follows. At each given moment most of the volume of the universe where 
the field <fi takes some particular value appear close to the centers of infloids created by the 
nonperturbative jumps by n(<f>)H/2iv. The new jumps occur each time H~ x independently of the 
previous history of the regions with a given 0. Therefore the leading contribution to the volume 
with be given by those rare centers of infloids where the field <fi jumps down by n(4>)H/2ir again 
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Figure 3: A schematic illustration which shows the number of infloids with given density and 
distribution of matter near their centers. 
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and again. That is why the typical volume weighted trajectories permanently go down with the 
speed exceeding the speed of classical rolling by n((p)H 2 /2ir, see eqs. (p8|) , (|4T|). 



One may visualize the resulting distribution of the scalar field in the following way. At some 
scale r the deviation of the field <ft from homogeneity can be approximately represented as a well 
of a radius r with the depth n(0)|^. Near the bottom of this well there is another well of a smaller 
radius e _1 r and approximately of the same depth n(0)|^. Near the center of this well there is 
another well of a radius e~ 2 r, etc. In particular, in the theory A0 4 /4 the depth of each well will 
be Mh^Ri._ Of course, this is just a discrete model. The shape of the smooth distribution of the 
scalar field is determined by the equation 



3H max H(f) 



dlnrH 2tt 




(49) 



which gives 



2 (r) « ^_ for r> H- 1 . (50) 

l-# max 2 (O)V?lrir# 

Note that <f>(r) ~ 0(0) for r < H^ 1 (there are no perturbations of the classical field on this scale). 

This distribution is slightly altered by the usual small perturbations of the scalar field. At a 
distance much greater than their wavelength from the center of the well these perturbations have 
the usual magnitude Thus, our results do not lead to considerable modifications of the usual 
density perturbations which lead to galaxy formation. However, the presence of the deep well 
(|50f) can significantly change the local geometry of the universe. 

In the inflationary scenario with V(4>) = ^0 4 fluctuations which presently have the scale 
comparable with the horizon radius ~ 10 28 cm have been formed at <p ~ 5 (in the units 
M p = 1). As we have mentioned already 3if max ~ 2y/6n ~ 8.68 for our choice of boundary 
conditions ||, and the typical nonperturbative jump down on the scale of the present horizon 
should be 3i7 max ~ 40 times greater than the standard jump, see eq. (P^). In the theory 
the standard jumps lead to density perturbations of the amplitude ^ ~ 2vjpi 03 ^ 5 . 
(in the normalization of 0). Thus, according to our analysis, the nonperturbative decrease of 
density on each length scale different from the previous one by the factor e should be about 
^ ~ -£f max ^y^0 4 ~ 2 ■ 10" 3 . This allows one to evaluate the shape of the resulting well in the 
density distribution as a function of the distance from its center. One can write the following 
equation for the scale dependence of density: 



1 dp 6V6Xk 4 

= ~H max — , (51) 



r r 
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where r is the distance from the center of the well. Note that = ^ (ln^)a in the theory 
|</> 4 ||. Here r corresponds to the smallest scale at which inflationary perturbations have been 
produced. This scale is model-dependent, but typically at present it is about 1 cm. This yields 



Ap = p(r) - p(r ) = 2ff max V6A r_ 

Pc p(r ) 57Ta/37t ~ r ' 

This gives the typical deviation of the density on the scale of the horizon (where In ^ ~ 60) from 
the density at the center: ^ ~ 750 • ^ ~ 4 • 10~ 2 . 

J Pc P 



It is very tempting to interpret this effect in such a way that the universe around us becomes 
locally open, with 1 — Q ~ 10 _1 . Indeed, our effect is very similar to the one discussed in 
p0| , ^T| , where it was shown that the universe becomes open if it is contained in the interior of a 
bubble created by the 0(4) symmetric tunneling. Our nonperturbative jumps look very similar 
to tunneling with the bubble formation. However, unlike in the case considered in pO , f2~l|j , our 
bubbles appear on all length scales. 



The results discussed above refer to the density distribution at the moment when the corre- 
sponding wavelengths were entering horizon. At the later stages gravitational instability should 
lead to growth of the corresponding density perturbations. Indeed, we know that density per- 
turbations on the galaxy scale have grown more than 10 4 times in the linear growth regime until 
they reached the amplitude ^ ~ 1, and then continued growing even further. The same can be 
expected in our case, but even in a more dramatic way since our "density perturbations" on all 
scales are much greater than the usual density perturbations which are responsible for galaxy 
formation. This would make the center of the well very deep; its density should be many orders 
of magnitude smaller than the density of the universe on the scale of horizon. This is not what 
we see around. 



This problem can be easily resolved. Indeed, our effect (but not the amplitude of the usual 
density perturbations) is proportional to H max , which is the maximal value of the Hubble constant 
compatible with inflation. If, for example, the maximal energy scale in quantum gravity or in 
string theory is given not by 10 19 GeV, but by 10 18 GeV, then the parameter i? max will decrease 
by a factor 10~ 2 . As we already mentioned, even smaller nonperturbative effects are expected 
in new inflation where -£f max is always many orders of magnitude smaller than 1. Inflationary 
Brans-Dicke cosmology in cases when the probability distribution P p is stationary also leads to 
negligibly small nonperturbative effects [13| . Thus it is easy to make our effect very small without 
disturbing the standard predictions of inflationary cosmology. However, it is quite possible that 
we will not have any difficulties even with very large n(4>) if we interpret our results more carefully. 
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8 Interpretation and possible improvements of the prob- 
ability measure 



An implicit hypothesis behind our interpretation is that we are typical, and therefore we live and 
make observations in those parts of the universe where most other people do. One may argue that 
the total number of observers which can live in domains with given properties (e.g. in domains 
with a given density) should be proportional to the total volume of these domains at a given 
time. However, our existence is determined not only by the local density of the universe but by 
the possibility for life to evolve for about 5 billion years on a planet of our type in a vicinity of a 
star of the type of the Sun. If, for example, we have density 10 -29 g • cm -3 in a small vicinity of 
the center of the infloid, and density 10 -27 g • cm -3 on the horizon scale, then the age of our part 
of the universe (or, to be more accurate, the time after the end of inflation) will be determined 
not by the density near the center of the infloid, by the large-scale density lCT 27 g • cm -3 , and it 
will be only about one billion years. 

Moreover, any structures such as galaxies or clusters cannot be formed near the centers of 
the infloids since the density there is very small. Indeed, on each particular scale the jump down 
completely overwhelms the amplitude of usual density perturbations. The bubble cannot contain 
any galaxies at the distance from the center comparable with the galaxy scale, it cannot contain 
any clusters at the distance comparable with the size of a cluster, etc. In other words, the center 
would be devoid of any structures necessary for the existence of our life. 

Thus, the naive idea that the number of observers is proportional to volume does not work at 
the distances from the centers which are smaller than the present size of the horizon. Even though 
at any given moment of time most of the volume of the universe at the density 10 _29 g ■ cm~ 3 is 
concentrated near the centers of infloids, the corresponding parts of the universe are too young and 
do not have any structures necessary for our existence. Volume alone does not mean much. We 
live on the surface of the Earth even though the volume of empty space around us is incomparably 
greater. 

One may argue, that the disparity between the age of the local part of the universe and its 
density appears only if one considers perturbations on a scale smaller than the horizon. Therefore 
it still may be true that we should live in the centers of huge bubbles, which have a shape fl52|) 
for r > Hq 1 , where Hq 1 is the size of the present horizon. If the cut-off occurs at r ^> Hq 1 , this 
may not lead to any observable consequences at all. However, if the cut-off occurs at r ~ Hq 1 , 
the resulting geometry may resemble an open universe with a scale-dependent effective parameter 

n(r) !. 

In order to make any definite conclusions about the preferable parts of the universe one should 
study probability distributions which include several other factors in addition to density. This 
should be a subject of a separate investigation. An additional ambiguity in interpretation of our 
results appears due to the dependence of the distribution P p on the choice of time parametrization. 
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Indeed, there are many different ways to define "time" in general relativity. If, for example, one 
measures time not by clock but by rulers and determines time by the degree of a local expansion 
of the universe, then in this "time" the rate of expansion of the universe does not depend on 
its density. As a result, our effect is absent in this time parametrization 0. The reason why 
the results depend on the time parametrization is deeply related to the properties of a self- 
reproducing universe. The total volume of all parts of such a universe diverges in the large time 
limit. Therefore when we are trying to find which parts of the universe have greater volume we 
are comparing infinities. There are some methods to regularize these infinities in a way that 
would make the final results only mildly dependent on the choice of time parametrization § . 
However, there are many such methods, and the final results are exponentially sensitive to the 
choice of the method ||. In this paper we used the standard time parametrization which is most 
closely related to our own nature (time measured by number of oscillations rather than by the 
distance to the nearby galaxies). But maybe we should use another time parametrization, see 
Appendix, or even integrate over all possible time parametrizations? Right now we still do not 
know what is the right way to go. We do not even know if it is right that we are typical and that 
we should live in domains of the greatest volume, see the discussion of this problem in M, B |i"5|. 



Therefore at present we would prefer to consider our results simply as a demonstration of 
nontrivial properties of the hypersurface of a given time in the fractal self-reproducing universe, 
without making any far-reaching conclusions concerning the structure of our own part of the 
universe. However, we must admit that we are amazed by the fact that the main fraction of 
volume of inflationary universe in a state with a given density p at any given moment of proper 
time t should be concentrated near the centers of deep spherically symmetric wells. We confirmed 
this result by four different methods, and we believe that it is correct. Until the interpretation 
problem is resolved, it will remain unclear whether our result is just a mathematical curiosity, or 
it can be considered as a real prediction of properties of our own part of the universe. At present 
we can neither prove nor disprove the last possibility, and this by itself is a very unexpected 
conclusion. Few years ago we would say that the possibility that we live in a local "center of 
the world" definitely contradicts basic principles of cosmology. Now we can only say that it is 
an open question to be studied both theoretically and experimentally. If somebody asks whether 
we should live in the center of the world, we will be unable to give a definite answer. But if 
observations show us that the answer is yes, we will know why. 
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Appendix 



Let us consider a different time parametrization, related to the proper time by local path depen- 
dent transformation: 



t-r(*) = f dsT(^(s)) , (53) 

where T(<p) is a positive function, and its argument in ([53]) is a solution of (|5]) with a particular 
realization of the white noise. The stochastic Langevin equation in this parametrization looks 
like follows: 



dr 3H(<t>)T((/>) 27^/2(0) w {0 > 

The branching diffusion equation in arbitrary time parametrization can be written as: 




^ 1 d ( H^) d ( H^((f)) 

dr pK9 ' } 2 d<f> \2irT 1 / 2 ((/))d(l) \27rT 1 / 2 (<f ) ) p 

+ Uml!k*™) +? $™ r) - (55) 



Its solution will generally be a stationary probability function with an overall constant expan- 
sion factor just like in ([14]) . The value of the constant Ai will depend on the parametrization. 

We can find the volume weighted slow roll trajectory of the inflaton field in arbitrary parametriza- 
tion very similarly to the approach used for proper time, but have to keep in mind that it is no 
longer true that Ai = df r H max . The result is: 



_ / v{4>) 

dr \\^H^)T{<j>) 



Ai 



y T{$)) 2vr 2 T(0) 



(56) 



Since the conventional (i.e. calculated under the comoving probability) amplitude of the quan- 
tum jumps generated during the typical time interval At ~ T(0)if _1 (0) in the given time 
parametrization is still given by the usual quantity H/2tt (see the Langevin equation ([54]) above) 
then the definition for amplification factor becomes: 
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n((f)) 



\ 



Ai 



2T{<j>) 2W{(t>) 
'T(0) ) H{<j>) 3HZ(<P) 



(57) 



In the particular case of the time parametrization T 
factor a(t) playing the role of time r, we get: 



H, which corresponds to the scale 



+ (Ai - 3) 



H 2 (<f>) 



2vr 2 



and 



n((f)) 



\ 



( W(<i>) \ 
{ 3# 3 (0) J 



+ 2 (Ai - 3) 



27^(0) 

3# 3 (0) 



(5* 



(59) 



Since Ai < 3, in this time parametrization the volume weighted slow roll (|58| ) is not faster 
but slightly slower than the conventional slow roll. As a result, most of the volume on the 
hypersurfaces of constant "time" r will be concentrated near the spherically symmetric hills 
(rather than wells) in energy density. However, the amplification factor is always very small. 

The change of time parametrization (|53D corresponds to one of the possible ways to choose 
regularization procedure for evaluation of divergent probabilities in eternally expanding universe 
0. Other types of regularization procedure were proposed in 0, In particular, the regular- 
ization scheme suggested in is essentially equivalent to choosing the T = H parametrization 
which we discussed above M. One can easily verify that in the limit <fi 0f r our equations for 
the T = H parametrization (0), (|59D yield the same results for the nonperturbative jumps as the 
ones obtained in [[?]]. As it is argued in ||, from the point of view of interpretation of our results 
it is not obvious that this regularization has any advantages as compared to a more intuitive 
and straightforward approach used in the main part of this paper. However, each regulariza- 
tion scheme and each time parametrization gives an additional interesting information about the 
structure of inflationary universe. Therefore we presented in this Appendix an extension of our 
results for the more general class of time parametrizations (0) . 
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